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We investigate the effect of induced massive radiation given off during the time of collapse of a 
massive spherically symmetric domain wall in the context of the functional Schrodinger formalism. 
Here we find that the introduction of mass suppresses the occupation number in the infrared regime 
of the induced radiation during the collapse. The suppression factor is found to be given by e~^ m , 
which is in agreement with the expected Planckian distribution of induced radiation. Thus a massive 
collapsing domain wall will radiate mostly (if not exclusively) massless scalar fields, making it 
difficult for the domain wall to shed any global quantum numbers and evaporate before the horizon 
is formed. 
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I. INTRODUCTION 

Black holes embody one of the greatest frontiers in 
theoretical physics - the intersection of quantum the- 
ory with general relativity. So far, attempts to confront 
this intersection and produce a single theory incorporat- 
ing both quantum mechanics and general relativity, have 
brought along with them several interesting properties 
(such as Hawking radiation, see solidified existing 
arguments (such as black hole entropy, see @]), and in- 
troduced unanswered paradoxes associated with the for- 
mation and evaporation of a black hole. These include 
the information loss due to the thermal nature of Hawk- 
ing evaporation (see 

One radical resolution of the information loss paradox 
was recently proposed in Q , where the collapse of a thin 
spherically symmetric domain wall (shell) was studied by 
means of the functional Schrodinger formalism (see also 
@, El)- Solving the Wheeler-de Witt equation for the 
collapsing shell, the authors of Q argued that (a) such 
collapse takes an infinite amount of time from the point 
of view of a distant observer, even when quantum fluctua- 
tions of the shell arc taken into account (but ignoring any 
back-reaction due to the emission of radiation), and (b) 
the shell emits pre-Hawking massless radiation prior to 
the formation of a trapped surface. This pre-Hawking ra- 
diation (PHR) could, they found, cause the shell to evap- 
orate completely in time, t cvap ~ Rs (l^r) (measured 
from the moment when the radiation begins to effectively 
carry away the mass of the shell). If (a) and (b) are both 
true, then no actual black hole can be formed in such a 
collapse, and, although the quantum state of PHR from 
the collapsing shell is entangled with the quantum state 
of the shell, information about the initial quantum state 
of the shell is continuously released during the process of 
evaporation. They found that the spectrum of PHR is 
approximately thermal [2l[ only at sufficiently large fre- 
quencies ujR s 3> 1, where R s is the Schwarzschild radius 
of the collapsing shell. The infrared part uiR s <C 1 of the 
spectrum, on the other hand, remains time-dependent 
and does not exhibit the Rayleigh- Jeans singular behav- 
ior characteristic of the Hawking radiation in finite time. 
The infrared part of the spectrum also depends on the 
initial quantum state of the shell and preserves all the 



information encoded in it. 

In the present paper we argue, to the contrary, that 
PHR will not prevent the formation of a trapped surface 
in the collapse of a spherical shell with mass M ^> Mp. 
First, we take into account quantum fluctuations of the 
Schwarzschild radius Rs ("quantum trembling" of the 
event horizon, see [l(|) of the shell, present the correction 
terms to the functional Schrodinger equation describing 
the quantum collapse of the spherical shell and show that 
such corrections make the collapse time finite and in fact 
rather short: t co n ~ Rs, counted from the moment when 
pre-Hawking evaporation starts being effective. Second, 
we show that the actual evaporation time i Gvap is signif- 
icantly longer than the time scale t CV ap ~ Rs (ijr) f° r 
pre-Hawking evaporation through emission of massless 
quanta, once one takes into account that the main con- 
tribution to the mass of the shell is given by constituent 
protons and neutrons. The latter implies that a typical 
collapsing shell has very large baryon/lepton quantum 
numbers, and in order for the shell to evaporate com- 
pletely prior to collapse, these quantum numbers must 
cither be erased by physical processes, or radiated away 
to infinity. Since global quantum numbers, such as B and 
L, are only carried by quanta of massive quantum fields, 
such as leptons and baryons, the relevant time scale to 
be compared to the total collapse time is the time of 
evaporation due to emission of massive PHR. 

We extend the analysis of Q to the case of PHR of 
massive particles and find the corresponding evapora- 
tion time. For masses of known leptons and baryons 
(including likely values of the neutrino masses) and as- 
trophysically interesting masses of the collapsing shell 
~M~ Hn~ 1 ^he total evaporation time is found 
to be significantly longer than the collapse time t co \\. 
We note how non-perturbative Standard Model baryon- 
numbcr violation softens, but does not eliminate, this 
conclusion. 

This paper is organized as follows. In the Sec. [TT] we 
construct the functional Schrodinger equation describing 
collapse of a quantum spherically symmetric shell and 
find corrections to this equation induced by fluctuations 
of the metric. We then compare classical and quantum 
collapse times for the shell. In Sec. IIIII we analyze mas- 
sive PHR emitted by the collapsing shell by means of 
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function Schrodingcr formalism, derive its spectrum and 
estimate its backreaction on the process of collapse. In 
Sec. IIVI we discuss the effect of baryon and lepton num- 
ber violation on the process of collapse. Finally, the Sec. 
IVlis devoted to the discussion and conclusions. 



II. ESTIMATES OF THE COLLAPSE 
TIMESCALE 

A. Collapse Time in Classical Physics 

Following let us consider a perfectly spherical do- 
main wall representing a spherical shell of collapsing mat- 
ter, so that the position of the shell is completely de- 
scribed by its radius, R(t). Our first goal is to check 
the well-known statement that from the point of view of 
a distant observer, the collapse takes infinite amount of 
time. Since we neglect any quantum effects, such as vac- 
uum fluctuations, the Schwarzschild radius of the shell, 
R s = 2GM . remains unchanged during the course to col- 
lapse, see [111 ] . 

Due to the spherical symmetry of the configuration 
and Birkhoff's Theorem, the metric can be naturally 
separated into two different regions. The metric of the 
space-time exterior to the collapsing shell is (see fTl| ) 
Schwarzschild: 



ds 2 



Rs 
r 



dt 2 



dr 2 



1 - R a /r 



where as usual 



dVt 2 = d6 2 + sin- 
In the interior of the shell, the metric is flat 
ds 2 = -dT 2 + dr 2 + r 2 d8 2 + r 2 sin 2 8d<j) 2 , r < R(t) (3) 

The interior time coordinate, T, is related to the asymp- 
totic observer time coordinate, t, via the proper time of 
an observer moving with the shell, r. The relations are 



r 2 dQ 2 , r > R{t) 
(1) 

(2) 



dT / ( dR 

^ = V 1+ [if 



and 



where 
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B(t) EE 1 



R(t.y 



(4) 



(5) 



(6) 



By taking the ratio of ((4]) and (|5|), the relationship 
between the interior time T and the asymptotic time t is 
found to be 



dT 
~dt 



VB + RJ 



B 



(7) 



where R T = dR/dr and R = dR/dt. 

Integrating the equations of motion for the shell gives 
an expression for the total mass of the shell [ll|, which 
is the integral of motion, 



M = -yi + R 2 + y/B + R 2 }4ir<jR 2 . (8) 

Here a is the surface tension of the shell. This expression 
for M is implicit since R s = 2GM occurs in B. Solving 
for M explicitly in terms of R T gives 



M = AnaR 2 [y/T+Rj - 2-kGgR\ . 



(9) 



When \R — R s \/R s <C 1, i.e. when the shell is close 
to its Schwarzschild radius, the velocity of the shell as 
observed by an asymptotic observer, is given by 

R ps —B . (10) 

Thus, to lowest order in (i?o — Rs), 

R(t) ~R S + (Ro - Rs)e~ t/Rs (11) 

where Rq = R(t = 0) is the initial position of the shell 
0- 

As we see, classical collapse of a spherically symmetric 
shell indeed takes infinite amount of time from the point 
of view of an asymptotic observer, although the distance 
between the shell and its Schwarzschild radius quickly 
becomes exponentially small. In what follows, we will 
count time starting from the moment where the behavior 
of the radius of the collapsing shell is approximately given 
by jnj|. 

B. Collapse Time in Quantum Physics 

We shall now consider how the classical estimate for 
the collapse time is modified once quantum effects (quan- 
tum fluctuations of the collapsing shell and of the metric) 
are taken into account. 

As shown in Q, the collapse time will remain infinite 
in the presence of quantum fluctuations in the position 
of the collapsing shell. Let us recall how exactly one 
comes to this conclusion. First, one assumes that the 
quantum dynamics of the system (the collapsing shell, 
the background space-time and the distant observer) are 
completely described by the Wheeler-de Witt equation 



if* = 



with the Hamiltonian 



if = if. 



sys 



if, 



obs i 



(12) 



(13) 



where ff b s is the Hamiltonian of the distant observer. 
Since the proper distance between the collapsing shell 
and the observer is very large, the quantum state of the 
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observer is not entangled with the state of the shell and 
the total wave function is separable: 



obs ' 



(14) 



One also expects that the state of observer \l> ob8 satisfies 
the usual Schrodinger equation 



.9* 



obs 



at 



where t is the observer's time coinciding with the usual 
Schwarzschild time. Therefore, the quantum state of the 
collapsing shell should also satisfy Schrddinger's equation 



at 



The Hamiltonian of the shell is Q 



H^ ys = BIl ■ MI + B{4iTfj,R Y 



(16) 



(17) 



where II = — i-^ is the canonical momentum, B is given 
by and /i = <r(l — 2irG<rRs) is a constant. In the near 
horizon limit R — > Rs, only the "kinetic" term survives 
in (|17p . and the "squared" functional Schrodinger equa- 
tion (fT6|) becomes especially easy to solve using tortoise 
coordinates t,u = R + Rs log \R/Rs — 1|: 



i:) 2 



if- 



dt 2 du 2 



* sys ss 0. 



(18) 



Its solution describes a wave packet propagating with 
the speed of light, with its width remaining fixed in the 
tortoise coordinates. Note however that the horizon is 
located at u = — oo, so that it takes infinite time for the 
wave packet to reach it. 

The considerations above do not, however, take all 
possible physical effects into account. For example, in 
deriving (|16|) . we assume a priori that neither the hori- 
zon scale Rs, nor the mass of the shell M = ^ fluctu- 
ates. However, both should fluctuate, at least when we 
try to measure their value with precision finer than the 
Planck scale, simply because the metric of space-time it- 
self strongly fluctuates at scales of the order of Planck 
length lp. Naively, it seems that these tiny Planckian 
fluctuations cannot be relevant for the physics discussed 
because the curvature of space-time at the scale R~Rs 
is significantly smaller than Mp. Recall however that 
the size of the shell R approaches the Schwarzschild ra- 
dius Rs exponentially quickly (sec (jTT|)) , so that the dif- 
ference R — Rs ~ i?oe~*/ i?s soon becomes comparable 
to the Planckian length lp. Therefore, at time scales 
t ^s> R s \og(Ro/lp) fluctuations of the metric must be 
taken into account. 

How can one then account for these fluctuations at the 
level of the Wheeler-de Witt equation in (|T2"j) ? By sepa- 
rating the metric fluctuations in the full Wheeler-de Witt 
equation from the fluctuations of the matter degrees of 



freedom associated with the collapsing shell and distant 
observer, one can write 



Hty = -16ttG- G ab , 



1 



cd " 



<JaW9cd 167TG 

+H sys + H ohs ) = 0, 

(19) 

where g a b is 3-metric of the space-time, R^ is 3- 



(15) curvature and G a 



b.cd 



9 



-1/2 



{g a c9bd 



gadgbc — 9ab9cd) 

is the local de Witt supermetric. The functional 
Schrodinger equation [l2| , as well as corrections to it [HI , 
can be constructed from (|19|) by looking for solutions of 
the form 



(20) 



with the "effective action" 5 expanded in powers of G: 

S = G- 1 S + S 1 + GS 2 + .... (21) 

The leading term So describes the classical background 
space-time and the classical collapsing shell, and equa- 
tions of motion derived by variating "the action" 5*0 are 
Einstein- Jacobi equations for the background. The term 
Si takes into account (a) the vacuum fluctuations of ge- 
ometry and (b) fluctuations of the shell and their backre- 
action on the background geometry. Variation of Si gives 
the functional Schrodinger equation discussed above. Fi- 
nally, the term S2 and higher order terms describe the 
backreaction of the vacuum fluctuations of geometry on 
the background as well as interaction between fluctua- 
tions of the shell and fluctuations of geometry. Equations 
of motion found by variation of S2 and higher order terms 
represent corrections to the functional Schrodinger equa- 
tion. 

The original Wheeler-de Witt equation (|T9)) is notori- 
ously hard to solve for many reasons, one of them be- 
ing that it requires regularization in the same way as 
does the Schrodinger equation describing behavior of the 
wave functional of relativistic quantum fields — behavior 
of a quantum field theory can be understood much more 
effectively by using second quantization and Fcynman di- 
agrammatics than first quantization and the Schrodinger 
equation. In fact, it is natural to expect that the very 
expansion in (|2ip breaks down at the vicinity of event 
horizon R — > Rs, when the behavior of the radius of the 
shell is approximately given by the expression in ([lip and 
where, as we shall see, the evaporation process seems to 
be the most effective. 

Fortunately, the relevant physics of the perturbation 
theory in powers of G can be understood using uncer- 
tainty relations. In particular, in order to estimate the 
collapse time one can notice that the horizon is expected 
to "tremble" at Planckian scales [l3 |. since the fluctua- 
tions of the metric g can be estimated as 



; 2 
' p 



(22) 
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where L is the length scale characterizing the curvature 



of the space-time |22| [14j. Correspondingly, the charac- 



teristic size of fluctuation of the Schwarzschild radius can 
be estimated as 



SR f 



/2 
l P 

Rs' 



(23) 



Recall again that at t ^> Rs the radius of the shell R is 
exponentially close to the Schwarzschild radius R s , and 
due to the uncertainty in (|23|) in determining the value of 
Rs, the observer at infinity will be unable to determine 
whether the actual event horizon is formed or not after 
the time 



t ~ R s log 



R s Rc 
' p 



(24) 



which can therefore be understood as the collapse time 
of the shell. 

Another way to reproduce the estimate of (j24|) is to 
use the fact that a distant observer can only determine 
the position of the collapsing shell by sending quanta to- 
wards it, having them scattered off the shell and measur- 
ing the corresponding cross sections (scattering and/or 
gravitational capture). Apparently, the maximal energy 
of the quantum, which the observer can send to the shell, 
is of the order of Mp. In the vicinity of the shell the 
blueshifted maximal energy of the quantum will there- 
fore be bounded by 



M P 



1 



R__ 

R 



(25) 



This energy is related to the uncertainty in determining 
the value of the position of the shell R(t) by 



5R- 



> 1, 



1 - 



showing that at 



t > 2R S log 



Ro 

77 



(26) 



(27) 



the observer will be unable to determine whether the 
shell has crossed its Schwarzschild radius. The estimate 
of |27p coincides with that of (|2"4")l , up to the logarithmic 
precision. 

As we see, once the fluctuations of metric are properly 
taken into account, the collapse time is no longer infinite 
(and in fact rather short taking into account how small is 
the value of the Schwarzschild radius Rs for astrophysi- 
cally interesting scales). 



III. PRE-HAWKING RADIATION OF A 
MASSIVE SCALAR FIELD 

In the process of collapsing toward its Schwarzschild 
radius, the shell emits PHR, losing mass. In principle, 



if this PHR has sufficient power, the shell could evapo- 
rate completely prior to formation of an event horizon, as 
argued in Q . To determine whether collapse can success- 
fully conclude, the collapse time scale as given in (|2"T|) has 
to be compared to the time scale of evaporation through 
PHR emission. Since the spectrum of the masslcss PHR 
is nearly thermal, with the temperature T being of the 
order of the Hawking temperature Q, the characteristic 
time scale for the evaporation of the collapsing shell of 
initial mass Mq through emission of massless quanta co- 
incides with the scale of Hawking evaporation for a black 
hole with Schwarzschild radius Rs = 2GMq 



MpR s = Rs 



Rs_\ 
h J 



(28) 



As long as Rs ^> lp, this time scale is significantly longer 
than the collapse time given in (|27p. and PHR cannot 
prevent the shell from collapse. 

Validity of this conclusion strongly depends on the es- 
timate in (|27|) for the collapse time. However, as we 
shall now argue, even if the estimate in (|27|) is ultimately 
wrong and an exact solution of the Wheeler-de Witt 
equation (fl~9| for a collapsing shell reveals a much longer 
(but finite) collapse time, the conclusion will hold un- 
der most (or perhaps all) astrophsically interesting con- 
ditions. The ultimate reason for its robustness is the fact 
that a collapsing shell made of protons, electrons and 
other constituents of the Standard Model carries global 
quantum numbers. These global quantum numbers are 
not known to be carried by massless particles. Since 
quanta of PHR with non-zero global quantum numbers 
have to be massive, the rate of their emission is expected 
to be exponentially suppressed by the usual Boltzmann 
factor [l| . Unless the Hawking temperature is at or above 
the mass of the particles, this Boltzmann factor will be 
a sizable suppression. 

In this Section, we shall expand the results of Q to the 
case of massive radiation, derive its spectrum and analyze 
its backreaction on the dynamics of a collapsing shell. 
In subsequent sections we shall examine the implication 
of the exponential suppression of massive modes on the 
likely evaporation pathways of astrophysical black holes. 



A. Functional Schrodinger Formalism for Massive 
Pre-Hawking Radiation 

Our goal in this section is to expand the results of Q 
to include the case of massive quanta radiated during the 
time of collapse. As a result, we will therefore carefully 
follow their development and notation. 

Let us consider a spherically symmetric configuration 
of a massive real scalar field $ which propagates in the 
background of the collapsing shell. The action for the 
scalar field is 



S 



d A x, 



-g\ {gTWdu* 



(29) 
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where g^ v is the background metric given by ([T} and (|3|) . 
The field <J> can always be expanded into a complete set 
of cigenmodes 



$ = 5>*(*)/*w 



(30) 



such that the Hamiltonian is a simple sum of terms. The 
total wavefunction then factorizes and can be found by 
solving a time-dependent Schrodinger equation of just 
one variable. 

From (JT|) and ([3]), one can split the action (|29p into 
interior and exterior parts 



S in =2TTjdtJ drr 2 T 
and 



y2 



■ + (<9r$r +m 2 1> 2 (31) 



S ou t^2TTldt I drr 

R(r) 



m 2 $ 2 



(32) 



where we used (JT]). It was shown in Q that ([7]) is given 
by 



dT 



R A 



(33) 



where /i = M/(4irfi) and /i = <r(l — 2iraGRs)- 

The regime of most interest is the one when the radius 
of the shell approaches the Schwarzschild radius R s . 

We see from (|3"3")l that in the near horizon limit, T ~ 
_B — > 0. Therefore the kinetic term for Si n diverges as 
(R— i?s) _1 i n this limit. The kinetic term in S ou t diverges 
logarithmically, so the Si n kinetic term is dominant term. 
Similarly the potential term in Si n vanishes while the po- 
tential term in S ou t becomes finite, so the potential term 
in S ou t dominates. Therefore we can write the action as 



S « 2tt / dt 



BJ 



drr 2 (d t <5>f 



drr 2 1 - 



R a 



{dr^f (34) 



m 2 I drr 2 <&<& 



where we have changed the limits of integration from R(t) 
to R s since this is the region of interest. 

Using the expansion in the modes (|30[) , we can rewrite 
the action as 



S 



dt 



- j^a k (t)M kk -d k >(t) + -a k (t)N kk *a k >(t) 



where d = da/dt, and M, N and P are matrices that are 
independent of R(t) and are given by 



M kk ,=4irl drr z f k {r)f k ,{r) 
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N fefe 



4-7T / drr 



(36) 

^r)/fc(r)&(r) (37) 
drr 2 f k (r)f k ,(r). (38) 



To simplify (1351) . we can write 



S-- 



dt\ 



--sa k (t)M kk ,d k ,(t) + a k (t)H kk ,a kl (t) 
rs 



where the matrix R is defined as 

RfcA- = N fcfc , + m 2 P 



(39) 



(40) 



From the action (f3"9"|) . we can find the Hamiltonian and, 
according to the standard quantization condition, 



n fc = -i 



d 



da k (t) 



(41) 



is the momentum operator conjugate to a k (t). From (|36[) 
and (|40p . we see that the matrices are Hermitian, there- 
fore it is possible to use the principle axis transformation 
to simultaneously diagonalize them (see Sec. 6-2 of [l5[ 
for example) . Then for a single eigenmode, the Hamilto- 
nian takes the form 



H — —- 



1 

2m 



db 2 



(42) 



where m and R denote eigenvalues of M and R, while b is 
the amplitude of the eigenmode. Hence the Hamiltonian 
for a single eigenmode takes (|42p the form of a harmonic 
oscillator with a time-dependent mass. 

To find the spectrum of the massive radiation, we must 
determine the occupation number of the quanta induced 
during the collapse. From the functional Schrodinger for- 
malism, the wave function ip(b,t) must satisfy 



dtp 



or from J42 



1 

2m 



1-5- 
R 



db 2 



-Rb 2 
2 



Re-writing (1431) in the standard form we obtain 



1 d 2 m ,., 
-2^W 2+ ^ )b 



ip{b,r]) 



. dipjb, rj) 
drj 



(44) 



where 



in , , , . 

— a k {t)P kk ia k >(t) 



(35) 



R 
mB 



B 



(45) 
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and 



77 = J dt'B. 



(46) 



Hence, instead of considering a time-dependent mass we 
consider the time dependent frequency. In (|45p we de- 
fined Wg = R/m. Here we will take the angular frequency 
wo to be 



w 



\J k 2 + m 2 



(47) 



which is the angular frequency for a relativistic particle. 
The exact solution to (j44"| is given in [lij 



V>(M) 





1/4 








im 






1 exp 


~2~ 









(48) 



where = dp/drj and p is given by the real solution of 
the non-linear auxiliary equation 



p m +UJ 2 (r])p = — 
P A 



with initial conditions 



p(0) = 



The phase a is given by 



^(0) = 0. 



1 z" 1 

a{r,) = -2j 



dij 

7Wr 



(49) 



(50) 



(51) 



Complete information about the radiation in the back- 
ground of the collapsing shell is contained in the wave- 
function (ggj) . 

If we consider that an observer at infinity will regis- 
ter quanta of the field $ at different frequencies, (|42|) 
tells us that the observer will interpret the wavefunction 
of a given mode at some later time in terms of simple 
harmonic oscillator states, {<p n }, at the final frequency, 
u>f. From (jllj) we see that in the near horizon, i.e. late 
time limit, B ~ e~ l l Rs . Therefore we can write the final 
frequency as 



u)f = woe 



t f /2R s 



(52) 



However, this is in terms of the conformal time, since the 
derivative in (|4"4")) is in respect to r\ not t. (|4T>|) tells us 
that the frequency in t is B times the frequency in 77, 
and at time tj, this implies that the observed physical 
frequency is then 



/*) = Buj 



f 



e-tf/X'Uf = w e-^ /2 ^ (53) 



where the superscript (t) on w refers to the fact that this 
frequency is with respect to t. The initial (t = 0) vacuum 
state for each modes is then simply the simple harmonic 
oscillator ground state 



1/4 



-majQb 1 2 



(54) 



where Wo = oj(t = 0). As mentioned previously, the 
number of quanta in each mode can be evaluated by de- 
composing (|48j) in terms of the simple harmonic oscilla- 
tor states and computing the corresponding occupation 
number per mode. The wavefunction for a given mode 
in terms of simple harmonic oscillator basis is given by 



where 



#,i) = ^c„(f)#) 



db<p*(b)ip(b,t) 



(55) 



(56) 



which is an overlap of the wavefunction at some later 
time t/j(b,t) with the simple harmonic oscillator ground 
state (|54[) . The occupation number at the final frequency 
uif is then given by 



(57) 



The occupation number in the mode b is then given by 
(see Appendix B of Q) 



N(t,e> f ) 



Pn 



(58) 



B. Numerical Results 

In the following plots, we have that the massless case 
(i.e. m = 0) is given by the blue line, m = 5R S is given 
by the purple line and m = 50Rs is given by the brown 
line. 

We have numerically evaluated the spectrum of mode 
occupation numbers at finite time and show the results 
in Figure [T] for different values of t/Rs- We see first that 
the occupation number of the massive radiation vanishes 
below = me"'/ 2Bs , as expected from (|4T|) . Above 
that, the spectra of the massless and massive cases are 
identical, since at this point the spectrum depends on 
uj^Rs only. It is also clear that the occupation number 
is non-thermal for low frequencies. In particular, as in Q, 
there is no singularity in TV at w = at finite time for 
the massless case. We also observe oscillations in N(oj) 
for both the massless and massive cases. 

We are interested in the effect that mass has on the 
spectrum of the pre-Hawking radiation given off during 
the time of collapse. We therefore plot \og(N) versus 
kR$ in Figure [5] for different fixed values of t/R s . Here 
we see that for small values of kRs the occupation num- 
ber is greater for smaller masses (largest for the massless 
case). However, for larger values of kRs the occupation 
of each of the masses begin to overlap, as expected from 
(f53|) . since as kRs increases w^ becomes dominated by 
k and the effect of the mass is then unimportant. The 
occupation number of the massive case is obviously sup- 
pressed when compared to the massless case. We wish 



t/Rs = 7, 8 

ln(N) 




F IG. 1: ln (JV) versus oj^Rs for t/R s = 7,8. Since o; (t) = 
Vfc 2 + 7n 2 e~ t//2i?s , for the massive case a/*' has a minimum 
value which is equal to me _t//2fls , as can be seen here. Note 
that the occupation for the massless case and the massive case 
overlap each other, since we are plotting versus w , which is 
the same for both cases after the minimum value is obtained. 
The occupation number grows as t/Rs increases 

to determine by how much the massive radiation is sup- 
pressed. 

To determine the amount of suppression, we compare 
the spectrum of mode numbers versus kRs with the oc- 
cupation numbers for the massless and massive Planck 
distribution, which are given by 

= ( 59 ) 

where uj = |fc| for the massless case and uj = y/k 2 + m? 
for the massive case, and where (3 is the inverse temper- 
ature. In Figures [2] and [3] we plot ln(A^) versus kRs for 
fixed values of t/Rs- Before we analyze the results we 
should note that the value of (3 in Figure [2] is different 
than that in Figure [3] This is because here we are plot- 
ting versus kR s not ui^Rs, hence we are ignoring the 
time dependence of the frequency, however keeping the 
time-dependence of the occupation number. Therefore 
each moment in time will have a different [3 associated 
with it. 

FigureOshows that for small values of kRs the numeri- 
cal distributions for the occupation number given in (|55|) 
are similar to the expected Planck distribution, shown 
here as dotted lines. However, as kRs increases, the two 
distributions differ from each other. This is due to the 
non-thermal nature of the radiation. 

Figure [3] shows that for small values of kRs the numer- 
ical distributions for the occupation number given in (|58[) 
are in agreement with the expected Planck distribution, 
again shown here as dotted lines. This is due to the fact 
that as t/Rs increases, the spectrum becomes more and 
more thermal for low frequencies, since we are ignoring 
the red-shift of the frequency. 

From (|59|) we see that for k « m the occupation 
number is exponentially suppressed by e~^ m . Due to the 
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t/R s = 10 

ln(N) 




FIG. 2: ln(JV) versus kR s for t/Rs = 10. The dotted lines 
correspond to the expected Planck distributions (|59[) . For 
small values of kRs the numerical values are similar, however 
as kRs increases the results differ. This is due to the non- 
thermal feature of the radiation for early times. 

t/Rs = 12 

ln(N) 




FIG. 3: ln(JV) versus kR s for t/R s = 13. The dotted lines 
correspond to the expected Planck distributions. Here the 
numerical results and expected Planck distributions are in 
agreement. This is due to the fact that as t/Rs increases the 
spectrum becomes more and more thermal. 



agreement in Figure [3] between the theoretical expected 
result (|59[) and the functional Schrodinger equation, we 
can also conclude that when considering massive radi- 
ation, the spectrum is exponentially suppressed by the 
mass. This result is in agreement with Hawking's predic- 
tion (see [H). 



C. Backreaction 

We arc finally ready to estimate the rate of pre- 
Hawking evaporation of the collapsing shell through the 
emission of quanta with mass m. The behavior of the 
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mass M of the shell is given by 



dM _ 1 dR s 
~dt ~ 2G~dT 



-AttR 2 



8ir 2 R s 



3/2 



(60) 

Taking into account (fTO]). with solution (JTTJ) in the 
near horizon regime (where the pre-Hawking emission is 
mostly effective), one can estimate 



time scale in (|62j) rather than that in ((28)) . For shells of 
solar mass I0 30 kg the time scale involved will be of the 
order of exp(10 19 ) years, apparently much longer than 
any physically interesting time scale in the problem. 



IV. EFFECT OF BARYON AND LEPTON 
NUMBER VIOLATION 



iy/8irmRs,o 



exp (AixmRs^ 



(61) 



at mRs.o S> 1. Here Rs,o is the Schwarzschild radius at 
the beginning of the collapse. 

As we explained in the Section lll B[ when the radius of 
the shell R becomes sufficiently close to Rs, so that one 
can use the expression in (JTTJ) , one naturally expects the 
perturbation theory in powers of G for the Wheeler-de 
Witt equation in (fT9]l to break down. This effect is ap- 
parently stronger than PHR, which is entirely contained 
in the order G° of the expansion corresponding to the ap- 
proximation of the functional Schrbdinger equation (note 
though that the backreaction of the PHR is not contained 
in the order G° ) . This effect is very well known in physics 
of black holes, where the effective Unruh temperature of 
the Hawking radiation increases when the observer ap- 
proaches the event horizon. We shall however assume 
that the conservative lower bound for evaporation time 
in (|61j) holds in all orders in G, since one expects higher 
orders of G to describe gravitational interaction (capture) 
of PHR with the shell as well as gravitational scattering 
of quanta of PHR on each other. 

Rewriting (|61[) in terms of the initial mass of the shell, 
M , 



tc 



M 2 

Rs< Mp 



exp 



/ AirmMa 



mMo 

We see that for M > Mp/Anm, 



/M ', 

iovap > Rs [jj- ) > Rs 



(62) 



(63) 



The condition M 3> Mp/iTrm, can be rewritten to place 
it more clearly into an astrophysical context: 



eV 

M » 10^ n M solar — 



(64) 



Note that among different global quantum numbers, 
baryon number B (and lepton number L) is of the most 
importance. One can naturally expect that the main 
contribution to the mass of the collapsing shell comes 
from masses of constituent protons and neutrons. As 
long as the pre-Hawking emission of particles carrying the 
baryon number remains suppressed, the shell cannot lose 
its mass effectively through the pre-Hawking emission of 
masslcss quanta. For such shells with masses larger than 
10 11 kg, the pre-Hawking evaporation time is given by the 



In the previous section, we have shown that the need 
to emit massive particles in order to radiate away global 
quantum numbers, such as the baryon or lepton numbers 
of the shell (or equally, gauge quantum numbers, such as 
the electro- magnetic charge of the shell), have the poten- 
tial to enormously extend the evaporation timescale of 
the shell. So far, however, we have not included the pos- 
sibility that physical processes which violate such global 
quantum numbers could allow a significant fraction of 
the mass of the shell to be radiated in masslcss quanta, 
leaving behind a less massive shell with a higher effective 
Hawking temperature. 

Indeed, the Standard Model contains a mechanism for 
changing the baryon and lepton numbers of the shell 
non-perturbative electroweak baryon and lepton number 
violation 

EMI- 

This process, which can, for example, 
convert 9 quarks into 3 anti-leptons [20l |. can allow the 
baryon number of our shell to relax. However, the process 
conserves B — L, indeed, it conserves B/3 — L e , — 
and B/3 — L T separately. (Here L, is the lepton number 
of the ith generation.) Thus, the minimum mass to which 
a shell of baryon number B (and lepton number L <C B) 
can evaporate by the emission of massless quanta is not 
~ Bitib, as one might have guessed, but rather 



B B 
M min (B) ~ —{m Vc + m„ + m„ T ) ~ — maxm,. 

3 3 i 



In other words 



M„ 



Mq maxj m Vi 
3 TUB 



(65) 



(66) 



We should note, that the electroweak non-perturbative 
baryon and lepton number violation is suppressed by the 



incredible e 87r/c 



o 3000 



at temperatures below the elec- 



troweak symmetry breaking scale. Thus in order for it 
to affect the evaporation of the shell materially, the elec- 
troweak symmetry must be restored. One might expect 
that this could happen in the region just above the sur- 
face of the shell as it reaches R s . We are going to address 
this issue in more details in the future work. 

If non-perturbative electroweak baryon number vio- 
lation does indeed convert the shell's baryons to anti- 
neutrinos, and thus reduce the mass by many orders of 
magnitude, neutrino flavor mixing could reduce it still 
further. Since v e , and v r are not mass eigenstates, it 
could in fact be that 



Mo mini m v - 



TUB 



(67) 
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where now i enumerates the (three or more) mass eigncs- 
tates. In this case we can rewrite 

(, 3mBM l ) 3 U Mm ^Z M °) ■ (68) 
evap b \(m™ n ) 2 M ) \ im B M 2 p J v ' 

In this case, i e vap 3> Rs only if 

^i^^M^^. (69) 

Interestingly, with m™ n ~ O.leV this is likely to mean 
that the largest observed black holes could decay by em- 
ploying these mechanisms. However, once again, it is not 
clear how neutrino oscillations would operate in this envi- 
ronment, and is probably enormously model dependent. 

We thus see that there is a potential for physics beyond 
the perturbativc standard model to act on the shell before 
it collapses. In this case, what would avert the collapse 
of the shell is not predominantly prc-Hawking radiation 
but the radiation of the high energy byproducts of the 
baryon-number violating interactions. 

V. CONCLUSIONS AND DISCUSSION 

As was recently demonstrated in Q , a spherically sym- 
metric collapsing shell generally loses energy due to pre- 
Hawking radiation. The spectrum of this radiation is 
nearly thermal for modes with momentum k 3> Rg 1 ~ 

^j- larger than the inverse Schwarzschild radius of the 
shell, while the occupation numbers of low momentum 
modes remain time-dependent until an event horizon is 
formed. If the time scale of this pre- Hawking evaporation 
is shorter than the collapse time, we expect this effect to 
prevent the formation of the event horizon and therefore 
lead to the resolution of the information loss paradox. 
The initial pure quantum state of the collapsing shell will 
evolve into a pure quantum state of the PHR, which con- 
tains the same amount information as the original state. 

The question whether such a scenario can be realized 
in practice is ultimately related to the hierarchy between 
two time scales in the problem: the collapse time for the 
shell and its evaporation time through emission of PHR. 
As we have argued, while classically the collapse time is 
infinite, taking quantum fluctuations of the space-time 
makes it finite and in fact rather small: 

icon w t + R s , (70) 

where to is the time necessary for the shell to collapse to 
the near horizon regime given in pip , irrelevant for the 
physics discussed, since prc-Hawking evaporation process 
is ineffective at t < to- 

The collapse time scale in ([70]) is to be compared to 
the total evaporation time for the collapsing shell due to 
emission of PHR. If one has a "realistic" collapsing shell 
made of baryons and lcptons, its mass is proportional to 



(approximately) the total number of baryons composing 
the shell and therefore to the total baryon number B of 
the shell. To evaporate the mass of the shell means to 
carry this baryon number away to infinity, and the rele- 
vant evaporation time scale should be associated to the 
process of emitting quanta of massive radiation. As we 
demonstrated, the spectrum of radiation is exponentially 
suppressed by the mass m of the emitted quanta, mean- 
ing that temperature of the black hole must be greater 
than that of the mass for that species to be radiated away. 
This would then imply that during most of the collapse 
of the massive domain wall, the radiation given off would 
be that of massless particles. For example, if a proton 
was to be radiated away during the collapse, the process 
would be suppressed by a factor of 

e -m e /T H _ e -10 10 _ 

In other words, it is difficult to have the collapsing matter 
evaporate before the trapped surface is created. 

Since massive radiation is only induced once the tem- 
perature of the black hole raises above that of the mass 
of the radiation, the first massive particle that would 
be radiated would be a neutrino, which carries lepton 
number, not baryon number. The lightest baryon, the 
proton, is at least ten orders of magnitude heavier, and 
thus its emission is horrendously suppressed. The shell 
can therefore lose its baryon number (and its mass) only 
if baryon number is violated, say by conversion into lep- 
tons. Such processes exist within the standard model, 
but are suppressed by approximately e -3000 , unless the 
electroweak symmetry is restored. Moreover, such pro- 
cesses would leave a third of the B— L to be carried by the 
heaviest neutrino, which may be too heavy to experience 
unsupprcssed PHR. Thus, within the Standard Model, 
effective radiation of the shell's mass may rely on all 
of electroweak symmetry restoration, non-pcrturbative 
baryon violation, neutrino oscillations, and PHR of neu- 
trinos happening sufficiently rapidly to evade horizon for- 
mation. Grand unified theories contain B-violating op- 
erators, and B (and even B — L) could be violated by 
gravitational operators above the Planck scale and short 
circuit some of this complication, but again, the rates 
of such processes in the extremely thin range of radius 
that is at GUT or Planck temperatures would need to 
be studied. Moreover, the radiation of the shell mass 
would proceed not so much by PHR as by the emission 
of high energy particles as a byproduct of baryon number 
violating interactions at high energy /temperature. 

Finally, let us note that the region where the pair pro- 
duction happens is located within vicinity of the shell 
(where the Schwarzschild metric describing the exterior 
of the shell is smoothly connected to the Minkowski 
space-time in the interior of the shell) , and while one par- 
ticle in the pair is escaped to infinity, the other particle 
from the pair enters the interior of the shell. Apparently, 
when the prc-Hawking evaporation process becomes ef- 
fective, the number density of particles in the interior of 
the shell grows, and eventually the original approxima- 
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tion of Minkowski space-time in <|3j> inside the shell breaks 
down. At this point, as one might naturally expect, (a) 
the space-time inside the shell will be described by the 
closed FRW metric and will collapse in finite time, and 
(b) density fluctuations in the gas of particles of PHR 
inside the shell will grow and collapse. This scenario ex- 
plains how exactly the formation of the trapped surface 
proceeds in the collapse of a spherically symmetric shell. 
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